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ABSTRACT 

Additional sources of CP violation in the MSSM may affect i?-meson mixings and decays, 
even in scenarios with minimal flavour violation (MFV). We formulate the maximally CP- 
^ \ violating and minimally flavour-violating (MCPMFV) variant of the MSSM, which has 
19 parameters, including 6 phases that violate CP. We then develop a manifestly flavour- 
covariant effective Lagrangian formalism for calculating Higgs-mediated FCNC observables 
in the MSSM at large tan (3, and analyze within the MCPMFV framework FCNC and other 
processes involving B mesons. We include a new class of dominant subleading contributions 
due to non-decoupling effects of the third-generation quarks. We present illustrative numer- 
ical results that include effects of the CP-odd MCPMFV parameters on Higgs and sparticle 
masses, the Bg and Bet mass differences, and on the decays Bg f^~^fJ'~, B^ tu and 
b s'-f. We use these results to derive illustrative constraints on the MCPMFV parame- 
ters imposed by DO, CDF, BELLE and BABAR measurements of B mesons, demonstrating 
how a potentially observable contribution to the CP asymmetry in the b ^ s'y decay may 
arise in the MSSM with MCPMFV. 

PACS numbers: 12.60.Jv, 13.20.He 
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1 Introduction 



Models incorporating supersymmetry (SUSY), such as the Minimal Supersymmetric Stan- 
dard Model (MSSM), contain many possible sources of flavour and CP violation. In partic- 
ular, the soft SUSY-breaking sector in general introduces many new sources of flavour and 
CP violation, giving rise to effects that may exceed the experimental limits by several orders 
of magnitude. The unitarity of the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing 
matrix suppresses flavour-changing- neutral currents (FCNC) and CP violation somewhat, 
thanks to the Glashow-Iliopoulos-Maiani (GIM) mechanism [1], to the extent that the 
soft SUSY-breaking scalar masses are universal. One possible solution to the flavour and 
CP problems is to ensure that the soft SUSY-breaking sector is fully protected by the 
GIM mechanism. This can be achieved within the so-called framework of minimal flavour 
violation (MFV), where all flavour and CP effects are mediated by the superpotential in- 
teractions corresponding to the ordinary Yukawa couplings of the Higgs bosons to quarks 
and leptons. In this framework, FCNC and CP-violating observables depend only on the 
fermion masses and their mixings, and hence the CKM mixing matrix V [2]. In such a 
scenario, all FCNC and CP violation observables would vanish in the MSSM if V were 
equal to the unit matrix 1. 

A minimal realization of MFV in the MSSM is obtained by assuming that all soft 
SUSY-breaking bilinear masses for the scalar particles, such as squarks, sleptons and Higgs 
bosons, are equal to a common value mo at the gauge couphng unification point Mgut, 
where Mgut might be the threshold for some underlying grand unified theory (GUT) 
based, e.g., on SU(5) or SO(IO). Likewise, the soft masses of the fermionic SUSY partners 
of the gauge fields, the gauginos, might also be equal to a common value mi/2 at Mqut 
and, in the same spirit, all soft trilinear Yukawa couplings of the Higgs bosons to squarks 
and sleptons could be real and equal to a universal parameter A times the corresponding 
Higgs-fermion-antifermion couplings. The Higgs supermultiplet mixing parameter fi and 
the corresponding soft SUSY-breaking term Bfi introduce two additional mass scales in the 
theory. However, minimization conditions on the Higgs potential can be used to eliminate 
these two last mass scales in favour of the electroweak scale Mz and tan/9 = Vu/vd, where 
Vu,d are the vacuum expectation values (VEVs) of the two Higgs doublets H^^ in the MSSM. 

It is well known that a minimal expansion of the above MFV framework is to allow 
the soft SUSY-breaking mass parameters mi/2 and A to be complex with CP-odd phases, 
thereby introducing two additional sources of CP violation in the theory. In this case, 
all FCNC observables, whether CP-conserving or not, still depend on the CKM mixing 
matrix V in such a way that they vanish if V is assumed to be diagonal, i.e., equal to 
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the unit matrix. However, the two new phases introduce the possibihty of CP violation 
in flavour-conserving processes even if V is real, and in general CP violation in FCNC 
processes may differ from CKM predictions. 

Here we go one step further, and ask the following question. What is the maximal 
number of additional CP-violating parameters and extra flavour-singlet mass scales that 
could be present in the MSSM, for which the above notion of MFV remains still valid, 
i.e., all FCNC effects vanish in the limit of a diagonal V? We call this scenario the max- 
imally CP-violating MSSM with minimal flavour violation, or in short, the MSSM with 
MCPMFV. As we will see in Section 2, there are a total of 19 parameters in the MSSM 
with MCPMFV, including 6 CP-violating phases and 13 real mass parameters. The pur- 
poses of this paper are to formulate the MSSM with MCPMFV, calculate the most relevant 
i?-meson observables, and explore the experimental constraints on the MCPMFV theoret- 
ical parameters, exploiting a manifestly fiavour-covariant effective Lagrangian formalism 
for calculating Higgs- mediated FCNC observables at large tan (3 that we develop here. 

At large values of tan/3, e.g. tan/5 > 40, one-loop threshold effects on Higgs-boson 
interactions to down-type quarks get enhanced [3-5], and so play an important role in 
FCNC processes, such as the K^-K^ mass difference, Bg-Bg and B^-Bd mixings, and the 
decays B — ^ Xg'j, B Kl^l^, Bg^ — > [6-15], and B ^ ru [16,17]. We present 

in this paper a manifestly fiavour-covariant effective Lagrangian formalism for calculating 
FCNC processes that follows the lines of the effective Lagrangian approach given in [12]. 
In addition, we include here the dominant subleading contributions to the one-loop Higgs- 
mediated FCNC interactions due to non-decoupling large Yukawa-coupling effects of the 
third-generation quarks. Based on this improved formalism, we compute FCNC observ- 
ables in constrained versions of the MSSM, where MFV has been imposed on the soft 
SUSY-breaking mass parameters as a boundary condition at the scale Mqut- We present 
numerical results for 5-meson observables in one example of the MCPMFV framework, 
from which illustrative constraints on the basic theoretical parameters are derived, after 
incorporating the recent experimental results from DO and CDF [18]. 

The paper is organized as follows: in Section 2, after briefly reviewing the MFV frame- 
work, we derive the maximal number of flavour- singlet mass parameters that can be present 
in the MSSM with MCPMFV at the GUT scale. All relevant one-loop RGBs are given in 
Appendix A. In Section 3, we present an effective Lagrangian formalism for Higgs-mediated 
FCNC interactions that respects flavour covariance. We also discuss the dominant sublead- 
ing effects at large tan /?, due to the large Yukawa couplings of the third generation. Useful 
relations which result from Ward identities (Wis) that involve the Z and VT-boson inter- 
actions to quarks are derived in Appendix B. Section 4 summarizes all relevant analytic 
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results pertinent to FCNC i?-meson observables. In Section 5 we exhibit numerical esti- 
mates and predictions for various FCNC processes, including the Bg-Bg and Bd-Bd mixings, 
and the decays Bg d A^"*^/^"; B Xg'-f, and B tu. We also illustrate the combined 
constraints on the theoretical parameters imposed by data from DO, CDF, BELLE and 
BABAR in one sample MCPMFV model. We summarize our conclusions in Section 6. 

2 Maximal CP and Minimal Flavour Violation 

In this section we derive the maximal number of CP-violating and real flavour- singlet mass 
parameters that can be present in the CP-violating MSSM and satisfy the property of MFV 
as described in the Introduction. 

The superpotential defining the flavour structure of the MSSM may be written as 

M^MSSM = U'^KQHu + D'^hdHdQ + E^XHdL + fiH^Hd , (2.1) 

where Hu,d are the two Higgs chiral superfields, and Q, L, f/*", and E'" are the left- 
and right-handed superfields related to up- and down-type quarks and charged leptons. 
The Yukawa couplings hu,d,e are 3x3 complex matrices describing the charged-lepton and 
quark masses and their mixings. The superpotential (2.1) contains one mass parameter, 
the /i parameter that mixes the Higgs supermultiplets, which has to be of the electroweak 
order for a natural realization of the Higgs mechanism. 

In an unconstrained version of the MSSM, there is a large number of different mass 
parameters present in the soft SUSY-breaking Lagrangian 

-jC^oft = ^ (MiSS + M2WW + M^rT + h.c.) + Q^MlQ + Z^MiZ + f/^M^f/ 
+ &yilb + E^mlE + Ml^HlH^ + MlHlHd + (^BfiH^Hd + h.c.) 
+ {u^Bi^QH^ + &s.dHdQ + E^s.,HdL + h.c.) . (2.2) 

Here Mi 2,3 are the soft SUSY-breaking masses associated with the U(l)y, SU(2)j;^ and 
SU(3)c gauginos, respectively. In addition, Mfj^ ^ and B^ are the soft masses related to 
the Higgs doublets and their bilinear mixing. Finally, '^'qldue 3x3 soft 

mass-squared matrices of squarks and sleptons, and g are the corresponding 3x3 soft 
Yukawa mass matrices ^. Hence, in addition to the yU term, the unconstrained CP- violating 
MSSM contains 109 real mass parameters. 

^Alternatively, the soft Yukawa mass matrices SLu.d.e may be defined by the relation: ia.u.d.e)ij = 
{i^u.d,e)ij {■^u.d,e)ij , where the parameters {A.u,d.e)ij are generically of order Msusy in gravity-mediated 
SUSY breaking models. In our paper, both definitions for the soft SUSY-breaking Yukawa couplings will 
be used, where convenient. 
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One frequently considers the constrained MSSM (CMSSM), which has a common gaug- 
ino mass mi/2, a common soft SUSY-breaking scalar mass mo and a common soft trilinear 
Yukawa coupling A for all squarks and sleptons at the GUT scale. The number of inde- 
pendent mass scales is greatly reduced since, even allowing for maximal CP violation, the 
free parameters are just mi/2, /i, mo, A and Bfi, where all but mo are complex variables. 
The phase arg fi may be removed by means of a global Peccei-Quinn (PQ) symmetry under 
which Hu and Hd have the same charges. Imposing the two CP-even tadpole conditions on 
the Higgs potential, one may replace fi = and Re (Bfi) by the Z-boson mass Mz and the 
ratio tan /3 = Vu/ f d of the VEVs of the Higgs doublets Hu^d, in the phase convention where 
Vu,d are real and positive. Linked to this, there is one extra CP-odd tadpole condition which 
can be used to eliminate Im (Bfi) in favour of maintaining the same phase convention for 
the VEVs, order by order in perturbation theory [19]. Thus, a convenient set of input mass 
parameters of the constrained CP-violating MSSM is 

tan/3(mt), mi/2(MGUT) , mo(MGUT) , ^(Mgut) , (2.3) 

where the relative sign of fi can always be absorbed into the phase definition of the complex 
parameters mi/2 and A. Thus, in addition to tan/3, this CP-violating CMSSM has just 5 
real mass parameters, two more than in its CP-conserving counter-part, namely the CP-odd 
parameters: Immi/2 and ImA. 

How can the general notion of MFV can be extended to this constrained CP-violating 
MSSM? In such a constrained model, the physical FCNC observables remain independent 
of details of the Yukawa texture chosen at the GUT scale. They depend only on the CKM 
mixing matrix V, the fermion masses, tan/5 and the 5 real mass parameters mentioned 
above. If the CKM matrix V were equal to the unit matrix 1, the FCNC observables 
would vanish, but flavour-conserving, CP- violating effects would still be present, associated 
with Immi/2 and ImA. Moreover, these parameters also contribute to CP- violating FCNC 
observables in the presence of non-trivial CKM mixing. Most noticeably, Immi/2 and ImA 
cannot generically mimic the effects of the usual CKM phase 6. 

We now consider how the above notion of MFV can be further extended within the 
more general CP-violating MSSM. To address this question, we first notice that under the 
unitary flavour rotations of the quark and lepton superfields, 

Q' = IJqQ, L' = IJlL, U'^ = VIjU^, D'^ = IJ*dD^, E'^ = IJ*eE^ , 

(2.4) 

the complete MSSM Lagrangian of the theory remains invariant provided the model pa- 
rameters are redefined as follows: 

h«,d —* ^l/^D ^u,d Uq , he he , 



5 



^Q,L,U,D,E ^Q,L,U,D,E^Q,L,U,D,E^Q,L,U,D,E , 

Siu,d a„,d Uq , He U^aeUL. (2.5) 

The remaining mass scales, /i, Mi^2,3, Mfj ^ and Bfi, do not transform under the uni- 
tary flavour rotations (2.4). In fact, it is apparent that the one-loop RGEs presented in 
Appendix A are invariant under the redefinitions in (2.5), provided the unitary flavour 
matrices Vq,l,u,d,e are taken to be independent of the RG scale. The effective Lagrangian 
formalism we describe in Section 3 respects manifestly the property of flavour covariance 
under the unitary transformations (2.4). 

It is apparent from (2.5) that the maximal set of flavour- singlet mass scales includes: 

Mi,2,3, M|^ ^, M^,L,C/,D,i? = ^W|,L,C/,Ai?l3, A„,,,, = I3 , (2.6) 

where the mass parameters fi and Bfi can be eliminated by virtue of a global PQ symmetry 
and by the CP-even and CP-odd minimization conditions on the Higgs potential. The 
scenario (2.6) has a total of 19 mass parameters that respect the general MFV property, 6 
of which are CP-odd, namely ImMi^2,3 and lmAu^d,e- 

We term this scenario the maximally CP-violating and minimally flavour-violating 
(MCPMFV) variant of the MSSM, or m short, the MSSM with MCPMFV. 

It is worth noting that, in addition to the flavour- singlet mass scales mentioned above, 
there may exist flavour non-singlet mass scales in the MSSM. For example, one could impose 
an unconventional boundary condition on the left-handed squark mass matrix Mg, such 
that 

MliMx) = M| Is + ml {h%) + ml (hlK) + mj {h%) (hlK) + . . . , (2.7) 

where Mx could be Mqut or some other scale. Evidently, there are in principle a consider- 
able number of extra mass parameters that could also be present in Mq(Mx), beyond 
the flavour-singlet mass scale Mq. In fact, these additional flavour non-singlet mass param- 
eters can be as many as 9 (including Mq), as determined by the dimensionality of the 
3x3 hermitian matrix Mq(Mx). The generalized boundary condition (2.7) on Mq(Mx) 
is in agreement with the notion of MFV for solving the flavour problem by suppressing 
the GIM-breaking effects, provided the hierarchy <C Mq is assumed. In particular, if 
these flavour-non-singlet mass parameters are induced by RG running, they may be 
generically much smaller than Mq. In this case, the will not all be independent of 
each other, e.g., in our MCPMFV scenario, the RG-induced flavour-non-singlet scales 
would be functionals of the 19 flavour-singlet mass parameters stated in (2.6). In general, 
a non-singlet mass parameter could either be introduced by hand or induced by RG run- 
ning of a theory beyond the MSSM with more flavour-singlet mass scales [20]. However, 
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since introducing -C Mq by hand has no strong theoretical motivation, we focus our 
attention here on the flavour-singlet MSSM framework embodied by the MCPMFV. 

Before calculating FCNC observables in the MSSM with MCPMFV, we first develop in 
the next section an effective Lagrangian approach to the computation of Higgs-mediated 
effects, which play an important role in our analysis. 



3 Effective Lagrangian Formalism 

Here we present a manifestly flavour- covariant effective Lagrangian formalism. This for- 
malism enables one to show the flavour-basis independence of FCNC observables in general 
soft SUSY-breaking scenarios of the MSSM. It will also be used in Section 4 to calculate 
FCNC processes in the MSSM with MCPMFV. 

To make contact between our notation and that used elsewhere in the literature [21], 
we redefine the Higgs doublets Hu^ as = $2 and Hd = iT2^l, where ti^2,3 are the 
usual Pauli matrices. We start our discussion by considering the effective Lagrangian that 
describes the tan /3-enhanced supersymmetric contributions to the down-type quark self- 
energies as shown in Fig. 1. The effective Lagrangian can be written in gauge-symmetric 
and flavour- covariant form as follows: 

= d'^nli^d^l + Ah,[$i,$2]) Q% + h.c, (3.1) 

where the superscript '0' indicates weak-eigenstate fields. In (3.1), the first term denotes 
the tree-level contribution and Aha is a 3 x 3 matrix which is a Coleman- Weinb erg- 
type [22] effective functional of the background Higgs doublets $12. We note that the 
one-loop effective functional Ah£;[$i,$2] has the same gauge and flavour transformation 
properties as hrf$J. Its analytic and flavour- covariant form may be calculated via 



(Ahrf), 



(27rj^ 



^.¥%lfl^— ^1 (3.2) 



1 
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where n = 4 — 2e is the usual number of analytically-continued dimensions in dimensional 
regularization (DR), 1^ stands for the x A^-dimensional unit matrix, Pl{r) = | [1~(+) 75] 
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Figure 1: Gauge- and flavour-invariant one-loop self-energy graphs for down-type quarks 
in the single-Higgs insertion approximation, with Hu = $2 cmd Hd = iT2^\. 



are the standard chirality-projection operators, and Cp = 4/3 is the quadratic Casimir 
invariant of QCD in the fundamental representation. The 8x8- and 12 x 12-dimensional 
matrices Mc and describe the squark and chargino-neutrahno mass spectrum in the 
background of non-vanishing Higgs doublets $1,2- 

It proves convenient to express the 8 x 8-dimensional chargino-neutralino mass matrix 
Mc in the Weyl basis (5, W^^'^'^, Hu, Hd), where Hu.d are SU(2)l doublets: Hu = 
{h^,h1)'^ and Hd = {h^,h'^)'^- In this weak basis, the Higgs-field-dependent chargino- 
neutralino mass matrix Mc[$i, $2] reads: 



Mc[$l,$2 





Mi 





^'^2 




\ 







M2I3 




-^,9^^{^r2) 








J- (b* 


O2 


H {iT2) 




\- 




^(7rf (zr2)$i 


-/i (2x2) 


02 


1 



(3.3) 



where g and g' are the SU(2)2, and U(l)y gauge couplings, respectively. Correspondingly, 
in the presence of non- vanishing Higgs doublets $1,2, the 12 x 12-dimensional squark mass 
matrix M^[$i, $2] is given by 



Qtf/ 



(M^)c7tQ (M^)f;t^ (M^)^t5 



(3.4) 



with 



= (Mj),,l2 + (hth,),,$i$t + (ht h„)., ($|$2 I2 - $2$ 
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(M2)^t^^^ = (M^,)., + (h„ht <|.t$2 + i 5,,. (^$1^1 - $t$2), 
(M2)5t5^, = (M^,),,. + (h,ht),,.<|.t$i - i5,,.^'2($t$i - 

(M2)^t5, = (M')~t^^ = (h„h;),,$fzr2$2, (3.5) 

where is the usual Kronecker symboL 

The form of the derived effective Lagrangian depends, to some extent, on the choice of 
renormahzation scheme. As usual, one may adopt the MS or DR schemes of renormaliza- 
tion. In general, the different schemes affect the holomorphic part of the Lagrangian at 
the one-loop level. Thanks to the non-renormalization theorems of SUSY, the Yukawa cou- 
plings h.u,d are not renormalized, and the wave functions of $1,2, QiL, um and dm remove the 
ultraviolet (UV) divergences of the one-loop corrections to the Yukawa couplings dm^iQjL 
and UiR^2QjL- The left-over UV-finite terms are not tan /3-enhanced and can be absorbed 
into the definition of hu^d, up to higher-order scheme-dependent corrections. Although the 
latter could be consistently included in our gauge-symmetric and flavour- covariant formal- 
ism, we ignore these small UV-finite holomorphic terms as they are higher-order effects 
beyond the one- loop approximation of our interest. 

By analogy, the gauge- and flavour- covariant effective Lagrangian for the up-type quark 
self-energies may be written down as follows: 

-/:;'ff[<l>i,$2] = iz°^(h„<l>n-2r2) + AK[^,,<!>2]),^Q% + h.c, (3.6) 

where Ah„[$i, $2] may be calculated from Feynman diagrams analogous to Fig. 1. As op- 
posed to the down-type quark self-energy case, these radiative corrections are not enhanced 
for large values of tan/3 and so are ignored in our numerical analysis in Section 5. 

The weak quark chiral states, m° ^ and (i^, are related to their respective mass eigen- 
states, Ul.r and dL,R, through the unitary transformations: 

< = U?ni, 4 = U^Vrfi, = Vlun, d^ = VUr , (3-7) 

where U^, U^''^ are 3x3 unitary matrices and V is the CKM mixing matrix. All these 
unitary matrices are determined by the simple mass renormahzation conditions: 

(^Cts[^i,<^2]) = -dRMddL + h.c, ^/:^ff[$i,<l>2]) = -urM^ul + h.c, (3.8) 
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where (. . .) denotes the value when the Higgs doublets $1,2 acquire their VEVs, and M^^rf 
are the physical diagonal mass matrices for the up- and down-type quarks. Imposing the 
conditions (3.8) yields [12] 



Q 



where 



^M,VtR, 

Vi 



R, 



1 + ^Uf 
1 + ^Uf 



h,iAh,[$i,$2] 



— M„R-i 

V2 



(3.9) 



U 



L ' 



U' 



Q 



(3.10) 



In (3.10) and in the following, the symbol 1 without a subscript will always denote the 
3x3 unit matrix. We observe that the unitary matrices U^, U^''^ can all be set to 1 
by virtue of the flavour transformations given in (2.4). The Yukawa couplings h^^^ are 
determined by the physical mass conditions (3.9). It is important to remark here [12] 
that these conditions form a coupled system of non-linear equations with respect to h^^d, 
since the Yukawa couplings also enter the right sides of (3.9) through the expressions Ild,u 
in (3.10). In addition, one should notice that the physical CKM mixing matrix V remains 
unitary throughout our effective Lagrangian approach. As we will see below and more 
explicitly in Appendix B, the unitarity of V throughout the renormalization process is a 
crucial property for maintaining the gauge symmetries through the Ward identities (Wis) 
in our effective Lagrangian formalism. 

We now consider the effective FCNC Lagrangian related to Higgs interactions to down- 
type quarks. From (3.1), we find that 



d,H 
efr 



- iai (1 + A 



ai 



Ul 



02 A 

h.c. 



02 



m2A 



0.2 



Yd, 



(3.11) 



where the individual components of the Higgs doublets $1,2 are given by 



$1, 



V2 



Vl,2 



X2 
01,2 



laio 



Moreover, the 3x3 matrices A^ ' , A^ ' and A^ ' are given by 



V2 



6 



501,2 



V2 



6 



5ai : 



(3.12) 



50t2 



(3.13) 
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where we have used the short-hand notation, = h^^ Ahc;[$i, $2], and suppressed the 
vanishing iso-doublet components on the LHS's of (3.13). In the CP- violating MSSM, 
the weak-state Higgs fields 01^2; 0,1,2 and 0]~2 are related to the neutral CP-mixed mass 
eigenstates Hi 2,3 [21,23], the charged Higgs boson H~ and the would-be Goldstone bosons 
and G~ , associated with the Z and W~ bosons, through: 

01 = On Hi, 02 = 02iHi, 

ai = cpG^ - spO-iiHi, a2 = s^G^ + cpOnHi, 

cpG- - s^H- , 02 = SfsG- + Cf^H- , (3.14) 



'1 



where = sin/3, = cos (3 and O is an orthogonal 3x3 Higgs-boson-mixing matrix. 

One may now exploit the properties of gauge- and fiavour-covariance of the effec- 
tive functional Ad[$i,$2] to obtain useful relations in the large-tan /3 limit. Specifically, 
Ad[$i, $2] should have the form: 

A,[$i,$2] = $Ifi + *Jf2 , (3.15) 

where fi^2 (^'^'I'^'i; '^*2'^2, '^'i'^'2, '^'2'^i j are calculable 3 x 3-dimensional functionals which 
transform as h^h^ or hj^h„ under the fiavour rotations (2.4). Given the form (3.15). it 
is then not difficult to show that in the infinite-tan /? limit [vi — > 0), 

limzv^(-^A,\ = ^(A,), lim(-^A,\ = ^(A,). (3.16) 

vi^O \ Oa2 I V2 '"1^0 \ 002 ' "^2 

Very similar relations may be derived for the up-type quark sector, but in the limit of 
vanishing tan/3. As we show in Appendix B, Ward identities (Wis) involving the W~ and 
Z-boson couplings to quarks give rise to the following exact relations: 

AT . .V2(^A.) ^ ^(A.). Ar . (^A.) ^ ^(A.>. (3.17) 

where v = \/vf + V2 is the VEV of the Higgs boson in the SM. Relations very analogous 
to those stated in (3.17) hold true for the up-type sector as well, i.e. A^" = A^^ = 
— -\/2 ( Au ) /v, where the extra minus sign comes from the opposite isospin of the up-type 
quarks with respect to the down-type quarks. 

For our phenomenological analysis in Section 4, we may conveniently express the general 
fiavour-changing (FC) effective Lagrangian for the interactions of the neutral and charged 
Higgs fields to the up- and down- type quarks u, d in the following form: 



C 



FC 



2M, 



w 
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(3.18) 



where the Higgs couphngs in the flavour basis = = = 1 are given by 
= ^VtR-(l + Aj^)v + ^VtR-A^V 

+ zOs. t, Vt R-i (l + A- - 1 A-) V , (3.19) 

g^,.-. = (giJJ^ (3.20) 
gU = —R-u'^t + ^R~'(l + At) 

Sf3 S/3 \ / 

+ z03.t^^R;^(l-A^^+t^A:^) , (3.21) 

Sluu = (gH,.n)^ (3.22) 

S^H-du = -t,VtR,i(l + AjQ + VtR-^Af , (3.23) 

g^-j. = -t,^Vt(l-(A^^)t)(R;i)t - vt(A^^)t(R;i)t, (3.24) 

and tp = tan (3. We note that the Higgs-boson vertex-correction matrices for the up-type 
quarks, A^^'^, A^^'^ and A^^'^ are deflned as in (3.13). 

The above general form of the effective Lagrangian Cpc extends the one derived in [12] 
in several aspects. First, it consistently includes all higher-order terms of the form 
(t^ mb/i/Mlugy)'^"^! which can become important in scenarios with large bottom-squark 
mixing [5]. Secondly, it does not suffer from the limitation that the soft SUSY-breaking scale 
should be much higher than the electroweak scale Mz- Speciflcally, SM electroweak cor- 
rections may be included in the Coleman- Weinb erg-type effective functionals Arf „[$i, $2], 
provided the theory is quantized in non-linear gauges [24] that preserve the Higgs-boson 
low-energy theorem (HLET) [25]. Finally, the effective Lagrangian Cpc implements prop- 
erly all the gauge symmetries through the Wis as discussed in Appendix B. 

The general FC effective Lagrangian (3.18) takes on the form presented in [12] in the 
single- Higgs-insertion approximation. In this case, the tan /3-enhanced threshold corrections 
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Figure 2: Two-Higgs-doublet model (2HDM) contribution to the one-loop self-energy graphs 
for down-type quarks in the single-Higgs-insertion approximation. 



«2 



, A^^, A^^ and (A^) are inter-related as follows: 

^^(A,) = A- = Aj^ = Af = (A:^ )■ , (3.25) 



V2 

where (A^) is given in the MSSM with MCPMFV by 



(3.26) 



and /(x, y, z) is the one-loop function: 



xy\xi{xy)+yz\n{yz)+xz\n{zx) 

I{x,y,z) = ^ . 3.27 

[x -y){y- z) [x - z) 

The ellipses in (3.26) denote the small contributions coming from the Feynman diagram in 
Fig. 1(c), which has the same flavour structure as the gluino-mediated graph in Fig. 1(a), 
i.e., this contribution is flavour- singlet in the single-Higgs-insertion approximation. We 
remark, finally, that in writing down (3.26) we have not considered the RG-running effects 
on the squark mass matrices between Mqut and MgusY- These effects are important, and 
are taken into account in our numerical analysis in Section 5. 

In addition to graphs involving SUSY particles, the two-Higgs-doublet model (2HDM) 
sector of the MSSM may also contribute significantly to the one-loop self-energy graphs 
of the down quarks. This contribution is shown in Fig. 2 and is formally enhanced at 
large tan since it is proportional to h^^. In the single-Higgs-insertion approximation, the 
2HDM contribution is given by 



V2 / 2HDM \ ^ HMu B*fi* 





+ 






+ 





(3.28) 
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This contribution turns out to be sub leading with respect to the Feynman diagram 1(b) and 
exhibits a very similar flavour structure. Beyond the single-Higgs-insertion approximation, 
the effective functional Ah^^°^[$i, $2] is calculated as 

where Mq[$i, $2] and M|^[$i, $2] are the 6x6- and 4 x 4-dimensional quark and Higgs- 
boson mass matrices in the background of non-zero $1,2- The 6 x 6-dimensional quark mass 
matrix is given by 

The Higgs-boson background mass matrix M^[$i,$2] receives appreciable radiative cor- 
rections beyond the tree level [19,21,23,26]. At the tree level, the 4 x 4-dimensional matrix 
M^[$i, $2] is given in the weak basis ($1, $2) by 

M5,|t.*d = i^^ff'^^ ^^ff-^A . (3.31) 



where 



(M|,)^t^^ = (M^);t^^ = (-5/. + ^$5<fi)l2 + (3.32) 

In the one-loop effective Lagrangian Cpc given in (3.18), the couplings of the Goldstone 
bosons and to quarks retain their tree-level form. This result is not accidental, but 
a consequence of the Goldstone theorem, which applies when the momenta of the external 
particles are all set to zero. However, the tree-level form of the Goldstone couplings gets 
modified when momentum-dependent (derivative) terms are considered. To leading order 
in a derivative expansion, one would have to consider the effective Lagrangian 



£9 = zQ^JZq^ + Ag'^)(<l>l(^$,) - (^$,)t$, 

+ Bg'^) ( (^$,)t - (^$,.) $t) ] + . . . , (3.33) 
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Figure 3: Dominant gauge- and flavour-invariant contribution leading to a modification of 
the tree-level Goldstone-boson couplings to quarks. 

where the dots denote analogous terms for the right-handed up- and down-type quarks Ur 
and dji. The first term depending on Zg is a functional of $1^2 for the left-handed quarks Ql- 
Such a term is not tan /3-enhanced and renormalization-scheme dependent. As mentioned 
above, these terms can be neglected to a good approximation. The effective functionals 
Aq"''*[$i, $2] and Bq-'''[$i, $2] are UV finite and include large Yukawa-coupling effects due 
to ht."^ In particular, this is the case for the effective functionals with i = j = 2. One 
typical graph of such a contribution is displayed in Fig. 3. Because of gauge invariance, 
analogous contributions will be present in the one-loop Z- and VT-boson couplings. All 
these effects are not enhanced by tan /3, and can be consistently neglected without spoiling 
the gauge symmetries of the effective Lagrangian £fc- 

In the next two sections, we present analytic and numerical results related to FCNC 
5-meson observables, using the effective Lagrangian (3.18) and including the 2HDM con- 
tribution (3.29). 

4 FCNC S-Meson Observables 

In this section, our interest will be in FCNC S-meson observables, such as the B^^-B^s 
mass differences AMb^^, and the decays Bs^d fJ'^P- , Bu tv and B Xg"^. 

^These effects tiave first been identified and studied in [27] witfiin the Standard Model. 
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4.1 AMb 



Our discussion and conventions here follow closely [12]. In the approximation of equal 
i?-meson lifetimes, the SM and SUSY contributions to AM^^ ^ may be written separately, 
as follows: 

AMs, = 2|(5j|iJ,^/=2|i?J)sM + (5J|ife^/=^|i?J)susY|, (4.1) 

where q = d,s and H^^'^ is the effective AB = 2 Hamiltonian. Neglecting the subdom- 
inant SM contribution, the SUSY contributions to the AB = 2 transition amplitudes are 
given by 

,1/2 



0.55 



(^.°|i.e^/-|i..W^2310ps-(|^)^(^ 

X "0.88(C2'^(°^) + c5^i^(2HDM)^ _ o.52(cf^(°P) + cf i^(DP))] , (4.2) 

where DP stands for the Higgs- mediated double-penguin effect. In addition, we have used 
the next-to- leading order QCD factors determined in [28-32], along with their hadronic 
matrix elements at the scale fi = 4.2 GeV: 

pLR ^ -0.58, P^^ = 0.88, Pf^^ = -0.52, P|^^ = -1.1. (4.3) 

The Wilson coefficients occurring in (4.2) are given by 

^SLL (DP) 

^1 



SRR(DP) 





3 

E 


„L _ „L _ 






167r^m^ 


3 

E 

i=l 




V2GfM^ 


Ml 


327r^mf,mg 


3 

E 

i=l 


„L _ „R _ 


V2GfM^ 


Ml 



^LR(DP) _ ■■"0''"q x ^ ^riioq^riioq ,n 
^^2 - - ,.9 / . ?72 ' l^-"^^ 

where the tan^ /5-enhanced couplings g^',^^ may be obtained from (3.18). Hence, the DP 
Wilson coefficients in (4.4) have a tan^ (3 dependence and, although two-loop suppressed, 
they become significant for large values of tan /3 > 40. 
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There are two relevant one-loop contributions to (i?°|if^^^^|i?°)susY at large tan/5: (i) 
the t-H^ box contribution to of the 2HDM type, and (ii) the one-loop chargino-stop 
box diagram contributing to Cf^^. To a good approximation, (^^r(2Hdm) given 
by |32] 



Cj-^'^^'^-^CWtan'/J. (4.5) 

In the kinematic region Mh± rrit, the above contribution can amount to as much as 
10% of the DP effects mentioned above. This estimate is obtained by noticing that the 
light-quark masses in (4.4) and (4.5) are running and are evaluated at the top-quark mass 
scale, i.e., ms{mt) ~ 90 MeV, md{mt) ^ 4 MeV [33]. The second contribution (ii) turns out 
to be non-negligible only for small values of the yU-parameter [32], i.e., for < 200 GeV. 

4.2 Bl^fi+fi- 

The leptonic decays of neutral B mesons, B^^ ^ /^^/^~; are enhanced at large values of 
tan/? [6-15]. Neglecting contributions proportional to the lighter quark masses m^,,; the 
relevant effective Hamiltonian for AB = 1 FCNC transitions is given by 

i/,^^=i = -2V2GFVt,V*(^CsOs + CpOp + Cio^io) , (4.6) 



where 



= mb (qPRb) (/i/i) , 

167r^ 

167r^ 

Oio = -^iqrPLb)m7,i,fi). (4.7) 

Using the resummed FCNC effective Lagrangian (3.18), the Wilson coefficients Cs and Cp 
in the region of large values of tan j3 are given by 



Cs 



27rm^ gH,qb9H,iif^ 



271171^ 1 A gH^gbdHr^f, 



where Cio = —4.221 denotes the leading SM contribution. In addition, the reduced scalar 

,P 

9k-,, = 9Hr,, = - tan/503.. (4.9) 



and pseudoscalar Higgs couplings to charged leptons g^^nn in (4.8) are given by 
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Here we neglect the non-holomorphic vertex effects on the leptonic sector since they are 
unobservably small. 

Taking into consideration the aforementioned approximations, the branching ratio for 
/i"*"/!" is found to be [8] 



(4.10) 



Ml 



Ml 



where q = d,s and TBg is the total lifetime of the Bq meson. Moreover, the form factors 

S,P,A 



Fl D A are given by 



'^MlFsg^^Cs,^ 
2 ^ mb + rriq 



n 



(4.11) 



Although the Wilson coefficient Ciq is subdominant for tan/? > 40, its effect has been 
included in our numerical estimates. 



4.3 B, 



TV 



There is an important tree- level charged-Higgs boson contribution to Bu tv decay [16, 
17]. It is not helicity suppressed and interferes destructively with the SM contribution [34]. 
The ratio of the branching ratio to the SM value is given by 



R 



Btu 



B{B- 



T U) 



1 + tan (3 ■ 



V 
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Mh± 



(4.12) 



where g^-j^ = — tan/? V"!" at tree level [cf. (3.23)], leading to the negative interference with 
the SM contribution. 



4.4 B Xs-f 

The relevant effective Hamiltonian for B — > Xg'y is given by 



h: 



b—*S'y 



eflf 



with 



4C I 

VtbV* { CMOM + c;(/i,)o;(/i,) + c',{fi,)o'M 



O2 



eiJif, 

16 7r2 
gs rrib 



16 vr 



sna^^F^'bL 



emi, 

16 7r2 
gsrrib _ 

TT^ Sr(^^luF'' bi 

16 TT^ 



(4.13) 



(4.14) 



1^ 



We closely follow the calculations of Refs. [35] for the branching ratio B{B — > ^s7) and 
the direct CP asymmetry in the decay. For the running c quark mass, we use rric (mP°i^) to 
capture a part of NNLO corrections [36]. We refer to, for example. Appendix B of Ref. [37] 
for the detailed expression of the branching ratio in terms of the Wilson coefficients which 
we are gong to present below. 

The LO charged-Higgs contribution is given by 



1 

3 



■>H'du' 



33 



V 



33 



V 



t 

23 



33 



V 



33 



V 



F^:iiy), (4.15) 



where y = fnf{Mw)/M'^±, the ratio of the top-quark running mass at the scale 
to the charged Higgs-boson pole mass. In the numerical analysis, we include the NLO 
contribution. Note that g^- = ~^/3^ and g]^- = —t/3 V' at tree level, see Eqs. (3.23) 
and (3.24). The functions Fyg*'*'^'' can be found in Ref. [37,38]. 
The chargino contributions are 

2 M?y 



(/^susy) 



i=l,2 



3 



(VtR-i)t3 Vt, + (VtR,i)t3 (Ci),2 {Cn)*,Mw ^(3), 



t \* 




23 

(M„R;:i)33 



72 



(Cl).2 {CRYaMw 



t \* 



-^7,8 \-^t. 



2 

m~ 



(4.16) 



where Xij = mf/m?-. We refer to [39] for the functions F^g* and to [40] for the chargino 
mixing matrices Cl,r and the stop mixing matrix f/*. 

Finally, the gluino contributions to the Wilson coefficients Cy^s are given by 



J2 



d\* 
L)i2 



i=l 



C|(/X; 



TTttc 



SUSY J 



V2Gf\M3\^X 



J2 ^^ iGi)*2 \ (Gi) 

Ms 



+ (<^^Ji3 



i3 



nib 



3/3(2;*) + 2/4(2;*) 



(4.17) 
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where A* = V33V23 = VtbV^* and Xi = |M3p/m~. The loop functions /i, 2,3,4(2^1) may be 
found in Ref. [41]. The Wilson coefficients for the primed operators O'^g can be obtained 
by the exchange L R and M3 —>■ M|: 



9V2GF\M,\^Xt 



X 



nib 



V2GF\M3\^Xt 



-yd \ 



3/1 (a;*) + -f2ixi) 



M* 
nib 



(4.18) 



In the above, C7 g^, the down- type squark-gluino- quark couplings G'l ^ are defined through 
the interaction Lagrangian (suppressing the colour indices) 



^dgd 



+d:[{Gi)iPn + {Gi)iPL] rf'dA, 



(4.19) 



where are the usual Gell-Mann matrices, i = 1, 2, . . . , 6 label the mass eigenstates of 
down-type squarks, and a = 1,2,3 the mass eigenstates of down- type quarks. The couphngs 
are given by the down-type squark mixing matrix as 



(G'i).« = u/'^n , (G^).« = -u/'^n . (4.20) 

The 6x6 unitary matrix diagonalizes the down-type squark mass matrix as 

^dt M|[/^~ = diag(m| , m| m| ), (4.21) 



where di is the lightest and d^ the heaviest. In the super-CKM basis, in which the down 
squarks are aligned with the down quarks and = = = 1, the 6x6 down-type 
squark mass matrix M~ takes on the form 



'-RR 



where the 3x3 submatrices are given by 



ML 



(ht h,) + C2P Ml 



-- + -s 



w 



(4.22) 
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ML = + |(h,hi) + CP Ml [-\sVj 1, (4.23) 

with h.d = ^ Mrf yl" R^^. As a byproduct of the chosen super-CKM basis, we observe the 
absence of flavour mixing in M~, for all h^-dependent terms, when R^^ oc 1. 

5 Numerical Examples 

For our numerical estimates of FCNC observables at large tan (3, we take the GUT scale 
to be the same as in the usual CMSSM with MFV, and a dedicated program has been 
developed to calculate the RG evolution from the GUT scale to the low-energy SUSY scale 
in the MCPMFV framework of the MSSM. For the Higgs mass spectrum and the mixing 
matrix Oai at the Msusy scale, the code CPsuperH [40] has been used. In the calculation of 
the flavour- changing effective couplings, only the leading contributions have been kept in 
the single-Higgs-insertion approximation, neglecting the EW corrections and the generically 
small flavour-off- diagonal elements of the squark mass matrices. 

In order to study the effects of CP-violating phases in the MCPMFV framework, we 
consider a CP-violating variant of a typical CMSSM scenario: 

|Mi,2,3| = 250 GeV, 

Ml^=Ml^ = Ml = Mlj = Ml = Ml = Ml = {im GeV)% 

= = = 100 GeV, (5.1) 

at the GUT scale with tan/3(MgusY) = 10, which corresponds to tan/? (m^"^*^) ^ 10.2. 
As for the CP- violating phases, we adopt the convention that $^ = 0°, and we vary the 
following three phases: 

$12 = $1 = $2 ; $3 ; = = = $Ae, (5.2) 

where, for simplicity, common phases $12 and are taken for the phases of Mi^2(^gut) 

and v4„ e(^GUT)5 respectively. We note that the phases of the gaugino mass parameters, 
$1 2,3, and the /i parameter, are unchanged by the RG evolution, whilst the phases 
of the elements of the matrix A„ ^ g could be significantly different at low scales from the 
values given at the GUT scale. This scenario becomes the SPSla point [42] when $1,2,3 = 0° 
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and $^UT ^ ^gQo_ Yiave found that MgusY varies between 530 GeV and 540 GeV, and 
^gut/10^^ GeV between 1.825 and 1.838 depending on the values of the CP-violating 
phases. 

We do not consider in this section the electric dipole moment constraints [43] on the 
MCPMFV parameter space of the MSSM. A systematic implementation of these constraints 
and their impact on the FCNC observables will be given in a forthcoming communication. 

5.1 Phases and Masses 

We first consider the (3,3) elements = (a/)33/(h/)33 at Msusy with / = u,d,e and 
fs = t, b, T. We find that the complex quantity Af^ can be written in terms of the complex 
Af and Mj at the GUT scale as: 

^/3(Msusy) ~ Cfj A;(Mgut) - M,(Mgut), (5.3) 

where the real coefficients C^J and Cj^' are functions of the Yukawa and gauge couplings. 
This expression is similar to that found in Ref. [44]. In general, cf^'"^ are much smaller than 
C^^'. Indeed, they are even smaller than ^^'^ with Cf"" < Cf'^. For A^, C^'= is not so much 
smaller than C^'-\ whilst Cf^ is negligible. This is because the strong coupling amplifies 
the infiuence of M3, while the large Yukawa couplings suppress those of the A terms via 
renormalization effects [44]. For the parameter set (5.1) with tan/5 = 10, we observe that 
the phases (Msusy) and (Msusy) are largely determined by $3, whereas the phase 
(Msusy) is niore affected by $1^2 than by This situation becomes different for 

larger values of tan/3, i.e. we find that C^^ becomes significant and C^"^ decreases when 
tan f] increases. 

In Fig. 4 we show sin<l>At, sin<l>yi^, and sin$yi^ for the parameter set (5.1) with 
tan /^(Msusy) = 10. In the left frames, we observe that $^1^ ^ and can be fully generated 
from $3 and $1,2, respectively, even when Au^d,e at the GUT scale are real, <I>^^"'" = 180°. 
Whilst the dependence of on $3 is negligible (solid line in the left-lower frame), the 
dependences of ^ on $12 can be sizeable (dashed fines in the left-upper and left-middle 
frames). In the right frames, the cases with $3 = 0° {^Att) "^12 = 0° i^A-r) are con- 
sidered, showing how large the A-term phases may become at the Msusy scale for real M3 
and/or real Mi and M2. When the gaugino masses are all real, | sin<I>^J and | sin$^J turn 
out to be 0.06 and 0.12, respectively, whereas |sin$^^| can be as large as 0.55. Some- 
what larger CP-violating phases are possible for and ^a^ when Mi and M2 are pure 
imaginary (see dashed and dash-dotted lines in the right-upper and right-middle frames of 
Fig. 4). Finally, there are no visible effects of $3 on 
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Figure 4: /n i/ie left frames, taking ^'^'^ = 180°, sin$^j (upper), sin^At (middle), and 
sin$^^ (lower) are shown as functions o/ $3 taking $12 = 0° (solid lines) and $12 taking 
$3 = 0° (dashed lines). In the right frames they are shown as functions 0/ taking 
$3 = 0° or $12 = 0°. For sin$^j anc? sin$^j^, three cases are shown: $12 = 270° (blue 
dash-dotted lines), 0° (black solid lines), and 90° (red dashed lines). For sin$^^, we set 
$3 = 0° as well. The parameters are taken as in Eq. (5.1) with tan/5(MsusY) = 10. 



We now discuss the effects of CP- violating phases on the masses of Higgs bosons, third- 
generation squarks and heavy neutrahnos and chargino. In the upper- left frame of Fig. 5, 
we show the absolute values of At^b,T as functions of a common phase $a/ = $1 = $2 = 
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Figure 5: The absolute values of At h ^ (upper-left) and the masses of the heavy Higgs bosons 
(upper right), sbottoms and stops (lower left), and charginos and neutralinos (lower right) 
as functions of a common phase $jv/ = $i = $2 = ^3- The solid lines are for = 180° 

and the dashed lines for = 0°. The parameters are listed in Eq. (5.1). 



for two values of 0° (dashed lines) and 180° (solid lines). In this case, one can show 

the absolute values squared depend only on the difference — ^^m'- 

\Af\^ ^ af-PfCO^{^T^ - ^m), (5.4) 

using Eq. (5.3), with aj ,Pf > 0. From Fig. 5, we observe that there is strong correlation 
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between |v4t^fc^T-| and the particle mass spectrum. This correlation is due to the phase- 
dependent terms Tr(ajjati) and Tr(a|jad) in dMfj^^^^/ dt and dM^ jjj^/dt. The fact that 
|M^^| decreases (increases) when Tr(a]^au) decreases (increases) explains the CP-odd phase 
dependence of heavier Higgs-boson masses, as can be seen from the upper-right frame of 
Fig. 5. The same correlation is observed for the heavy chargino and neutralinos in the 
lower-right frame of Fig. 5, since a decreased (increased) value of leads to smaller 

(larger) values of We find that the variations in the masses of the lightest Higgs boson 
Hi and the lightest neutralino x? amount to 2 GeV and 3 GeV, respectively. The CP-odd 
phase dependences of Mg, M^, and M|, at the scale Msusy can be understood similarly. 
Here the (3,3) components of the mass matrices decrease (increase) when Tr(aj^a„) increases 
(decreases). For the chosen value of tan /^(Msusy) = 10, the (3,3) component of shows 
the largest effect, since dM^/dt contains 2Tr(aJ^a„) compared to Tr(a|^a„) -|- Tr(a)^ad) in 
d'M.^/dt and 2Tr(a]^ad) in dWl\)/dt. Furthermore, we note that ii ~ Ir and hi ~ hi- From 
these observations, one can understand the qualitative CP-odd phase dependence of the 
stop and sbottom masses, as shown in the lower- left frame of Fig. 5. 



5.2 Effects on AM^^ and AM^, 

In the upper- left frame of Fig. 6, we show the SUSY contribution to AMb^ in units of 
ps~^ as a function of tan/3(MgusY) for three values of the common phase, namely <I>m = 0° 
(solid line), 90° (dashed line), and 180° (dash-dotted line). The horizontal line is for the 
measured value: AM^f^ = 17.77 ± 0.10 (stat.) ± 0.07 (syst.) ps"^ [18]. We observe that 
the SUSY contribution can be larger than the current observed value for = 180° when 
tan/3 is large. Indeed, for $Af = 180° (90°), we find tan/3 < 44 (48), whereas there is no 
restriction on tan/3 for = 0°. 

The SUSY contribution jg suppressed by ml /ml with respect to Q^^^^^^^i 

[see Eq. (4.4)]. The |C^^^°^^| is comparable to |Cf^'"''°^^ |, while the 2HDM contribution, 
^LR(2HDM)^ becomes less important as tan/3 increases. The dip of the coupling 
for $A/ = 180° (upper-right frame) at tan/3 ^ 45 is due to the fact that the three Higgs 
bosons become degenerate and cancel other contributions. Beyond this point, Mh^ ~ 
decreases rapidly while Mh^ ~ 110 GeV remains nearly unchanged. 

In the upper-left frame of Fig. 7, we show the SUSY contribution to AMb^ in units 
of ps~^ as a function of tan/?(MsusY)5 using the same line conventions as in Fig. 6. The 
horizontal line is for the measured value: AM^^ = 0.507 ± 0.005 ps~^ [45]. We observe 
that the SUSY contribution is always smaller than the measured value, although it does 
exhibit a strong dependence on the CP- violating phase $m- The dips at tan/? ^ 45 



25 



,— , 10 



10 



< 



10 



10 



EXP 



MEASURED 



20 



30 



0^=180" 



O =0 



40 



50 



10 

Q 

^ 10 

U 
— 10 

10 

10 

10 

10 



-2 



tanp 



-3 



-4 



-5 



-6 



-7 



-8 



10 



20 



30 



40 



50 



tanp 




-2 



-3 



-4 



-5 



-6 



-7 



-8 



10 



20 



30 



40 50 

tanp 



Figure 6: The SUSY contribution to AMb^ in units o/ps~^ (upper-left) and the relevant 
couplings in the other three frames, as functions o/ tan/3(MsusY); for three values of the 
common phase: = 0° (solid lines), 90° (dashed lines), and 180° (dash-dotted lines). 
We fix = 0° and the parameters are taken as in Eq. (5.1), except that here we choose 

Mi^E = 200 GeV so as to avoid a very light or tachyonic fi state for large tan/3. In the 
upper-left frame, we show the currently measured value as the horizontal line. 

($M = 180°) and tan/5 ~ 49 = 90°) arise for the same reason as in the AM^^ case. 
The dominant contribution comes from q^^^^^^^ ^ g^^d (^^^^(^p) jg suppressed by m^/m^. 
The value of ICg^''^^'*! is smaller than that of |Cf'"'"''^^'' |. Finally, as before, the 2HDM 
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Figure 7: The SUSY contribution to AMb^ in units o/ps"^ (upper-left), and the relevant 
couplings in the other three frames. The line conventions and the parameters are the same 
as in Fig. 6. 

contribution (^^^(^hdm) i^g^^j^gg ^ggg significant for large values of tan (3. 
5.3 Effects on Bs fi^fi^ 

In the upper- left frame of Fig. 8, we show the branching ratio B{Bs l^^fJ' ) as a function 
of tan/5(MsusY) using the same line conventions as in Fig. 6 for three values of the common 
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Figure 8: The branching ratio B{Bs Ai^/W") in the upper-left frame and the relevant 
couplings in the other three frames, in units of GeV~^ as functions o/tan/3(MgusY)- The 
line conventions and the parameters chosen are the same as in Fig. 6, except that the two 
horizontal lines in the upper-left frame are for the SM prediction and the current upper 
limit at 90 % C.L. 



phase ^M- = 0° (solid line), 90° (dashed line), and 180° (dash-dotted hne). The two 
horizontal lines in the upper-left frame are for the SM prediction and the current upper 
limit at 90 % C.L., namely 7.5 x 10~® [18]. We observe that the branching ratio changes 
substantially as $a/ varies. Specifically, for $m = 180° (90°) 0°, we find that the present 
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upper limit on B{Bs fj^^fi ) imposes the upper limit tan/5 < 34 (38) 42. 

The phase dependence of the branching ratio comes from that of the couplings Cs and 
Cp [see (4.8)], which are shown in the upper-right and the lower-left frames, respectively. 
We find that \Cs\ ^ \Cp\, since On ~ 0«i ~ and Mh^ ~ Mh, [cf. (4.8) and (4.9)]. 
We note that, for $a/ = 180°, B{Bs — > /i''"/i~) can be smaller than the SM prediction for 
tan/3 < 24. This is because the Higgs-mediated contribution Cp cancels the SM one Cio, 
as shown in the lower-right frame of Fig. 8, in which the factor mfe/(mfc + m^) [cf. (4.11)] 
has been suppressed in the label of the |/-axis. 

5.4 Effects on tu 

The recent BELLE and BABAR results for the branching ratio B[B^ r~u) are [46,47] 

B{B- ^ r-p)™ = (l.79^°i^ (stat)+Sf (syst)) x 10-\ (5.5) 
B{B- ^-p)BABAR ^ (^-^ 2 ± 0.4 (stat) ± 0.3 (bkg syst) ± 0.2 (other syst)) x 10^^ , 

which lead to B{B- r-u)^^^ = (1.4±0.43) x 10"^ Combining the BELLE and BABAR 
results with the SM value B{B- r^uf^ = (1.41 ± 0.33) x 10"^ obtained by the global 
fit without using B{B~ t~u) as an input [48], we have the following 1 a range for the 
ratio to the SM prediction 

R^^l = 1.0 ± 0.38. (5.6) 

In the upper-left frame of Fig. 9, we show possible values of this ratio in the MSSM 
with MCPMFV, together with the experimental range given in (5.6), as functions of tan jS 
for three representative values of the common phase $m and for = 0. The three 

thin arrows at the bottom indicate the positions where the ratio vanishes at the tree 
level without including threshold corrections for = 180°, 90°, and 0° (from left to 
right). Beyond the minimum point, the charged Higgs-boson contribution dominates over 
the SM one. It rapidly grows as tan^ /3 initially and then goes over to tan^ /3 due to the 
threshold corrections. For each displayed value of $m, we find two regions of tan/5 where 
the experimental value of B{B'^ t~u) is obtained. One region is at tan/5 < 25 (27) 29 
for $jv/ = 180° (90°) 0°, and corresponds to the case where the charged Higgs-boson 
contribution is a small 'correction' to the SM term. The second region is at tan/5 ~ 
41 (46) 48, for $a/ = 180° (90°) 0°, and corresponds to the case where the charged Higgs- 
boson contribution dominates over the SM term. We note that the locations of these 
second allowed regions would not be estimated correctly if the threshold corrections were 

•^This range is different from that used in [49] due to the new BABAR resuh [47]. 
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Figure 9: The ratio Rbtu (upper-left), the charged-Higgs boson mass in GeV (upper-right) , 
and the real (lower-left) and imaginary (lower-right) parts of the coupling (g^t^j^)i3/Vi3 = 
(g^* j^)3i/\4fe as functions of tan f] for three or four values of taking ^^^"^ = 0°. 
The experimentally allowed 1-a region is bounded with two horizontal lines in the upper- 
left frame. The straight line with a tag 'Tree' in the lower-left frame shows the tree-level 
coupling. The parameters are the same as in Fig. 6. 

not included. These regions are actually excluded by the Bg — > constraint discussed 

previously. 

The tree-level vanishing points are also indicated in the upper-right frame as intersec- 
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tions of the Mh± and tan j3 x M^± lines. We observe that the resummed threshold effects 
enhance the charged Higgs-boson contribution when $m = 180° and suppress it when 
$M = 0°. As can be seen from the lower-left frame of Fig. 9, for $a/ = 90°, the tan/5- 
dependence of Rbtu becomes rather similar to the tree- level one. However, as displayed in 
the lower-right frame of Fig. 9, there is a non- vanishing contribution from the imaginary 
part of the coupling (gJ^tj„)i3/Vi3. 

5.5 Effects on B ^ Xs^ 

The current experimental bound on B{B X^'-f) with a photon energy cut of > E^ut = 

1.6 GeV is [50] 

B{B Xsif^^ = (3.55 ± 0.24+°:°^ ± 0.03) x 10"^ (5.7) 

Our estimate of the SM prediction based on the NLO calculation is 3.35 x 10~^, which is 
about 1 a larger than the NNLO result, (3.15 ± 0.23) x 10""^ [36]. In Fig. 10 we show the 
branching ratio B{B ^ ^sl) and the direct CP asymmetry Aqp{B — s> Xs'~f) as functions 
of tau/S. In the upper-left frame, we include only the charged-Higgs contribution, which 
increases the branching ratio. The larger contribution in the high-tan /3 region is due to 
the decrease of the charged Higgs-boson mass. In the upper- right frame of Fig. 10, we 
add the contribution from the chargino-mediated loops. This contribution largely cancels 
the charged-Higgs contribution, when < 90°. Instead, if $m is larger than ~ 90°, 
the chargino contribution interferes constructively with the SM one, resulting in a rapid 
increase of the branching ratio as tan (3 grows. This behaviour can be understood from the 
fact that the dominant contribution to Cy § comes from the last term of Eq. (4.16), which is 
proportional to ~ e**^t/c/3, and the branching ratio is proportional to its real part, namely 
cos$At/c/3. We recall that the phase at the low-energy scale can largely be induced 
by non- vanishing even when vanishes (see the upper frames of Fig. 4). In the 

lower-left frame of Fig. 10, we show the full result including the contribution of the gluino- 
mediated loops, which is non-vanishing in the presence of flavour mixing in the down-type 
squark mass matrix. We find that it is numerically negligible for the parameters chosen. In 
the same frame, as well as in the upper-right one, we show the case of the common phase 
= 60°, in which there is a nearly exact cancellation between the chargino and charged- 
Higgs contributions, and all the tan/3 region considered is compatible with the current 
experimental bound. This observation is also apparent in the left panel of Fig. 11. In the 
lower-right frame of Fig. 10, we show the direct CP asymmetry for several combinations of 
(<|)GUT ^ <|)^^)^ finding that it can be as large as ~ —4 %, when <I>m = 60°. 
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Figure 10: The branching ratio B{B Xg'-f) as a function o/tan f3 for several values of the 
common phase $a/ = $i = $2 = ^3 ond The region allowed experimentally at the 

2-a level is bounded by two horizontal lines. In the upper-left frame, only the charged-Higgs 
contribution is added to the SM prediction. In the upper-right and lower-left frames, the 
SUSY contributions are included. The direct CP asymmetry Aqp{B Xg'-f) is also shown 
in the lower-right frame for several combinations of ($5^"'" , $Af)- The parameters are the 
same as in Fig. 6. 

To illustrate the strong dependences of the branching ratio and the CP asymmetry on 
the common phase $M; we show them as functions of <I>a/ for four values of tan/5 in Fig. 11. 
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Figure 11: The branching ratio B[B ^ ^sl) (^^ft) and the CP asymmetry Aqp{B Xg'-f) 
(right) as functions of for four values of tan jS taking = 0°. The region allowed 

experimentally at the 2-o level is bounded by two horizontal lines in the left frame. In the 
right frame, points satisfying this constraint are denoted by open squares. The parameters 
are the same as in Fig. 6. 

The region allowed experimentally at the 2-a level is bounded by two horizontal lines in 
the left frame. In the right frame, points within this region are denoted with open squares. 
We observe that the branching ratio is quite insensitive to tan/5 around = 60°, whereas 
the CP asymmetry can be as large as ±5 % for points within the current 2-a bound on the 
branching ratio. For comparison, we note that the experimental range currently allowed is 
0.4 ±3.7 % [50], implying that the new contribution in the MSSM with MCPMFV could be 
comparable to the present experimental error, and much larger than the SM contribution, 
which is expected to be below 1%. Finally, it is important to remark that, in the absence of 
any cancellation mechanism [43] , EDM constraints severely restrict the soft CP-odd phases 
in constrained models of low-scale SUSY, such as the constrained MSSM. In a forthcoming 
paper, however, we will demonstrate in detail, how these constraints can be considerably 
relaxed in the MSSM with MCPMFV. 
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6 Conclusions 



In this paper we have formulated the maximally CP-violating version of the MSSM with 
minimal flavour violation, the MSSM with MCPMFV, showing that it has 19 parame- 
ters, including 6 additional CP-violating phases beyond the CKM phase in the SM. As 
preparation for our discussion of i?-meson observables, we have developed a manifestly 
flavour- covariant effective Lagrangian formalism, including a new class of dominant sub- 
leading contributions due to non- decoupling effects of the third-generation quarks. We have 
presented analytical results for a range of different i?-meson observables, including the Bg 
and Bd mass differences, and the decays Bg f^^fJ'", B^ tu and b — » 57. We have 
presented numerical results for these observables in one specific MCPMFV scenario. This 
serves to demonstrate that the experimental constraints on 5-meson mixings and their de- 
cays impose constraints, e.g., on tan/5, that depend strongly on the CP-violating phases in 
the MCPMFV model, most notably on the soft gluino-mass phase in the specific example 
studied. 

In summary, on the one hand, our paper introduces a new class of MSSM models of 
potential phenomenological interest and develops an appropriate formalism for analyzing 
them, and on the other, it presents exploratory numerical studies of the constraints im- 
posed by experimental limits on S-meson observables. In view of the large number of the 
theoretical parameters in the MSSM with MCPMFV, we leave for future work a more com- 
plete exploration of its parameter space, including the correlation with other experimental 
constraints, e.g. those imposed by limits on electric dipole moments. 
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A Renormalization Group Equations 



Here we list all relevant one-loop renormalization group equations (RGEs) for the gauge 
and Yukawa couplings [51], as well as for the soft SUSY-breaking mass parameters of the 
general MSSM [52,53]. Defining the RG evolution parameter t = In(QV^GUT)) 

we may 

write down the one-loop RGEs as follows:^ 
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"^Our results are in agreement with [53]. 
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+ 2M^>ht + 2aeat + ^ glTi (YM^) l^j , (A.17) 

where gi is the GUT-normahzed gauge couphng, which is related to the U(1)y gauge 
couphng g' of the SM through gi = a/5/3 g'. In addition, the expression 

Tr{YM^) = Mjj^ - Ml^ + Tr (^Mj - - 2M^ + + M|) (A.18) 

is the Fayet-Ihopoulos D-term contribution to the one-loop RGEs. It can be shown that 
(iTr (yM^)/(it oc Tr(yM^), i.e., the expression Tr(yM^) is multiplicatively renormal- 
izable. As usual, the GUT scale is determined by the boundary condition: 5'i(Mgut) = 
fl'2(^GUT) = fi'sl^GUx)- We note, finally, that the one-loop RGEs listed above are invariant 
under the unitary fiavour transformations given in (2.5). 



B Z- and VF^-Boson Ward Identities 

In the absence of gauge quantum corrections, the Z- and boson couplings to quarks 
obey the following tree- level Wis [54]: 



rf^'(g,P,P-g) + r^''^^'(g,P,P-g) = (B.i) 

iV^^^^^q, p.p-q) + 2r^^"'^(g, p,p-q) = (B.2) 



LW 



w 



(p) Pl - Vud' Pr ^dd'ip - q) 



where = a/1 — is the cosine of the weak mixing angle and T"*-"^"* = \{ — \) and 
Qu(d) = I (~|) the weak isospin and electric charge quantum numbers for the u and 
d quarks. In (B.I) and (B.2), T.ffr{p) are quark self-energies describing the fermionic 
transition /' — > /, with f = u,d and /' = u',d'. In addition, T^^f'{q,p,p — q) and 
r|f ^"'^(g,p,p — q) are vertex functions that describe the interaction of the Z- and W~^- 
boson to quarks, respectively. The momenta q^ of the gauge bosons are defined as fiowing 
into the vertex, while the momentum fiow of the quarks follows the fermion arrow, where 
always denotes the outgoing momentum. 

In general, virtual strong and electroweak gauge corrections to the Z- and ly^-boson 
vertices usually distort these identities, through terms that depend on the gauge-fixing 
parameter, e.g, ^. One possible framework in which these identities can be enforced is 
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the pinch technique [55], leading to analytic results that are independent of C,. Recently, 
this approach has been extended to super Yang-Mills theories [56]. We ignore the gauge 
quantum corrections in our phenomenological analysis, since they are rather small. 

In the limit 0, the Wis (B.l) and (B.2) simplify considerably. Let us first consider 
the WI involving the Z boson and its associate would-be Goldstone boson G^. Since the 
vertex function T^^^' (q,p,p — q) has no IR singularities in the limit — > 0, the WI (B.l) 
takes on the much simpler form 



^9n 



Mzc^ 



Tl 



(B.3) 



Decomposing the quark self-energies Ilffi{p) with respect to their spinorial structure, 

Effip) = Sf^,(p^)^Pi + J:^{p')^Pn + ^ff,ip')PL + ^f,}ip')PR- (B.4) 



we may rewrite (B.3) as follows: 

r«"//'(o,p,p) = 



T/ 



w 



^if,{p')PL - PR^f*f{p'' 



(B.5) 



Considering the proper normalizations determined by the relations given in (3.13), it is 
possible to make the following identifications in the effective potential limit 0: 

Sf_,,(0) = Ujh;(A;)U{, P,r^°//' (0,0,0) = --^Ujh^Af U{P,, (B.6) 

where the unitary matrices U^'^ take care of the weak to the mass basis transformations 
as given in (3.7), with = and = V. Then, the simplified WI (B.5) implies 
that 



V2 



(B.7) 



which is the relation assumed in Section 3 [cf. (3.17)]. 

We now turn our attention to the WI involving the boson and the associated 
would-be Goldstone boson G"^. In the effective potential limit q^,p^^ 0, we obtain 



G+ud 



(0,0,0) 



w 



V.,,S^„,(0)Pi - V„,,S^*(0)P^^ 



:b.8) 



Employing the definitions (3.13) and taking the weak-to-mass basis rotations of the quark 
states into account, we find the relations: 



P,,T 



G+ud 



(0,0,0) 



Pz.r^"""(o,o,o) = u^Th,A:^"u^Pi 



:b.9) 
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From the simplified WI (B.8) and its Hermitean conjugate, we then derive that 

Af = -^i^u), AT = ^{A,), (B.IO) 

which is in agreement with (3.17) and the discussion given below. We note that the unitarity 
of the radiatively-corrected CKM matrix V lies at the heart of deriving the relations (B.7) 
and (B.IO). 



C CPsuperH Interface 



To solve the RGEs given in Appendix A, we have considered the following input parameters: 
• The gauge couplings at the scale Mz: 



5 g'\Mz 



a2{Mz 



(C.l) 



3 Air ' ' 47r 

where g{Mz) = e{Mz)/sw and g'{Mz) = e{Mz)/cw with a,^{Mz) = e\Mz)/in. 
The evolutions of ai^2 from Mz to m^°^'^ are determined by [57] 
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On the other hand, a3(m^°^'^) has been obtained by solving the following equation 
iteratively [58] 



) = aj'(Mz)-6oln 
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where = -(11 - 2Nf/3)/2tt, hi = -(51 - 19 NF/3)/A7r^, and 62 
5033 Nf/9 + 325 N'^/27) /GAir^ with Np = 5. 
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The masses of the quarks and the charged leptons at the top-quark pole-mass 
scale m^"^'^. In particular, the top-quark running mass at m^"^'^ is obtained from: 

The CKM matrix V is assumed to be given 



mt(m^°^'^) = 171^°^^ / 1 + 4:as{m^°^^) /3n 

t the same scale 
are given by 



at the same scale m^"^'^. Then, in general, the complex 3x3 Yukawa matrices at 
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(C.4) 
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in the flavour basis = = = I3. The diagonal quark and charged-lepton 
mass matrices are given by 





= diag 


rriu 






1 pole\ 
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Given 01,2,3(^^4'°^'^) ^u,d,e{^t°^'^), the evolution from m^"^'^ to the scale Msusy 
have been obtained by solving the SM RGEs. Here the SUSY scale Msusy has been 
determined by solving 



iteratively, where m? = max(m^ + mf, im? 



=A^SUSY 
^3 



max[m?((5^) , fnf^iQ'' 



(C.6) 



U3 



m?) and m? = maxim 



ml). For m|^ ^^(Q^), we have taken the (3, 3) component of the corresponding 

mass matrix as 

"'Q3,C/3,£'3,L3,£;3 

At the scale Msusy, the Yukawa matrices match as 



(3,3) 



(C.7) 



hd,e(Ms+jsY) 



h«(MgusY)/sin/5(MsusY), 

hrf,e(MgusY)/cOS/3(MsuSY) 



(C.^ 



and, finally, the evolution from Msusy to Mqut have been obtained by solving the 
MSSM RGEs. 

The 19 flavour- singlet mass scales of the MSSM with MCPMFV, which are parame- 
terized as follows: 



U/f |o«*l,2,3 14 |p**A 

M'^i,2,3| e ' ' , \^u,d,e\ e 



u,d,e 



,U,D,L,E ; 



M; 



(C.9) 



These are inputed at the GUT scale Mqut, which is defined as the scale where the 
couplings gi and g2 meet. Any difference between 5'3(Mgut) and 5'i(Mgut) may be 
attributed to some unknown threshold effect at the GUT scale. 



By solving the RGEs from the GUT scale Mqut to the SUSY scale Msusy, we obtain: 

— Three complex gaugino masses, |Mj| e**'((5 = Msusy)- 

— Three 3x3 complex Yukawa coupling matrices, hu,d,e{Q = ^susy)- 
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Three 3x3 complex a-term matrices, SLu,d,e{Q = ^susy)- 
The soft Higgs masses, M^^^ jj^iQ = Mgusv)- 

The complex 3x3 sfermion mass matrices, ^q^u.D.L.EiQ = ^susy)- 



The inputs for the code CPsuperH are: 



t^nPimf'), Mir, /i(MsusY), Mi,2,3(Msusy) 



"^03,03,1)3,^3,^3 (^susy), A(Msusy), A(Msusy), A,(Msusy) • (C.IO) 
The ratio of the vacuum expectation values at mJl°^^ is related to that at Msusy by [21] 
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The gaugino mass parameters are directly read from the results of the RG running, the 
sfermion masses are given by 

"^03,C/3,i'3,i3,^3(^SUSY) 



Mq,c/,d,l,£;(^susy) 



(3,3) 



(C.13) 



and the A parameters, including their CP-violating phases, by 

[a/ (Msusy)] (3,3) 



Af{M, 



SUSYj 



[h/(MsusY)] 



(C.14) 



(3,3) 



The yU parameter and charged Higgs-boson pole mass can be obtained from 

M^^ (Msusy) and M^^ (Msusy) by imposing the two CP-even tadpole conditions, T,^^ = 
T^2 = [21]. The tadpoles can be cast into the form 



^</.i(02) = ^^1(2) /i?(2) + ^'2(1) Remi2 + t'l(2) 

where 

^1(2) 



^1(2) f?{2) + ^(^3 + A4)t'2{l) 



m 



pole 2 



^1(2)^1(2) (C.15) 



(C.16) 



The quantities fif ^ and Aj are given by 



/^?,2 
X; 



A. + A«(mn + Af (mr^^). 



(C.17) 
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where 



2(1)/ polo\ 
2(1)/ polc\ 



167r2 
3 



M? 



pole 2 



m 



pole 2 



M? 



pole 2 



m 



m 



pole 2 



(C.18) 



The couphngs Aj, A^^^(m^°^°) and Aj-^''(mf°'*') may be found in Ref. [21]. The squared 
absolute value can be determined from {T^^/v2 — T^^/vi) = 0, which does not depend 



on Rem^2! since 



tb 



(C.19) 



with 



X, 



X, 



167r2 
3 



l/iJ'IAI^n. 





pole 2 

ml 



M? 



m 



pole 2 



167r2 



|/i*l'ln- 



pole 2 



M? 



rri- 



pole2 



(C.20) 



We note that the phase of the /i parameter, is not renormalized. 



Once l/ip is found, Remi2 can be obtained from T^^ = or T^^ = 0- With Rem 



12 



known, the charged Higgs-boson pole mass can be obtained by solving the following equa- 
tion iteratively: 



spcp 2 



(C.21) 



For the explicit form of Hh+h-^ we refer to Ref. [59]. We note that, for large tan/3, 
Rem^2/"5/3C/3 — Mjj^ — Mfj^ — M| at the tree level. Finally, after imposing the CP-odd 
tadpole condition Im(i?/i) = 0, we use = Rem^j calculate the 2HDM contribu- 
tion (3.28), by noting H^Ha = -$1$2- 
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